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Abstract. This paper describes the analysis used to determine the neutral current Tt'^ production in 
MiniBooNE in bins of momentum. Additionally, a measurement of the relative coherent production 
of tt'^s is discussed. The coherent production rate is found to be (19.5 ±1.1 (stat) ±2.5 (sys))% of 
the total exclusive neutral current tt" production rate. 

INTRODUCTION 

Neutral current K^— production is a potential major background to the Vg appearance 
signal that MiniBooNE is looking for in its test of LSND [1]. As such the main objective 
of this analysis is to measure the rate of TZ^ production so that misidentification in the 
Vg oscillation sample can be determined as a function of reconstructed Vg energy in 
the charged current quasi-elastic mode. As a result the initial product of the NC Tt^ 
analysis effort is not an absolute cross section. Instead it is a measurement of the total 
TZ^ production in bins of momentum, and a measurement of the coherent production 
which effectively fixes the angular distribution. The analysis described here was used 
in the recent oscillation analysis reported by MiniBooNE [2*] to fix the Tt^ production 
in the Monte Carlo based on the observed rates in data. In addition, the dynamics 
of neutrino induced, neutral current n'^ production is of interest in its own right. The 
analysis shown here builds upon the work previously shown at Nulnt04 [[30. 

DETERMINING tt^ PRODUCTION IN MOMENTUM 

The event selection begins with a set of pre-cuts that exactly match the pre-cuts used 
in the electron neutrino selection. The event must have only a primary event without 
evidence of a secondary event consistent with a muon decay (or Michel electron). This 
eliminates the vast majority of charged current interactions. The event must have 
more than 200 hits in the main tank. This is well above the Michel endpoint. The 
event must have fewer than 6 veto hits. This eliminates more that 99.9% of all cosmic 
rays. Additionally, all events must be in the 1.6 jUs beam spill window, although by the 
time the all other cuts are applied this is essentially all that remains. The determined 
production rates are all relative to these pre-cuts. Therefore, if one is interested in 
computing a cross section from these numbers, it would be important to understand 
the inefficiency of these cuts (for example from the overlap of two neutrino interactions. 
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or of a single neutrino interaction with a cosmic ray) and the effective target volume, 
which is largely set by the veto cut. 

The analysis cuts are based on the reconstruction which fits each event with muon, 
electron and hypotheses. Each fit produces a likelihood, and the log of the ratio 
of different likelihood hypotheses are used for particle identification. In the first stage, 
electron-//fce events are selected by applying the cut log(^^/^e) > 0.05. Next the n'^- 
like events are selected with log(^;i:/=^e) < 0. When coupled with a reconstructed 77 
mass cut, about the mass, a very clean sample of events (signal to noise ratio ~30) 
is selected with an overall efficiency of about 40%. 

The n*^ candidate events are divided into bins of reconstructed n'^ momentum, and 
the Monte Carlo (MC) is used to unsmear the data, correcting for momentum smearing 
and inefficiency. A matrix is formed by dividing MC events into bins of true momentum 
verse reconstructed momentum and counting events over background in each bin. Each 
matrix element is divided by the total number of events in all reconstructed bins with 
the same true momentum range. This matrix, which is well conditioned and largely 
diagonal, is inverted to form the unsmearing matrix. The data vector is formed by 
dividing the candidate events into the same reconstructed bins and subtracting the 
background in each bin according to the signal to noise ratio for the corresponding 
MC reconstructed range. The unsmeared, or corrected, data rates are the product of 
the unsmearing matrix and the data vector. Figure [U shows an absolutely normalized 
comparison of the raw MC prediction to the corrected data distribution. The ratio of 
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FIGURE 1. Results of the unsmearing in bins of momentum. The blue points show the corrected 
piO momentum distribution and the red points show the raw Monte Carlo momentum distribution. This 
comparison is shown absolutely normalized. 
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FIGURE 2. The reweighting function, which is used to correct the tt" rate of the Monte Carlo in 
bins of momentum. This function is just the ratio of the distributions shown in Figure [1] 



these two distribution forms a correction function (Figure |2l) wliicli is used to reweight 
events as a function of true momentum, in the MC. 

By construction, the reweighting fixes the discrepancy between data and MC in re- 
constructed momentum. Additionally, it also improves agreement in many of the key 
kinematic distributions. Figure [3] shows the relatively normalized data to MC compari- 
son for both the raw and corrected MC. The kinematic distributions shown are the cosine 
of the 77 opening angle, the 7 energies and the 77 mass. The opening angle and energy 
comparisons show marked improvement, while agreement in the mass distribution is 
largely unaffected by the reweighting. 



RESONANT AND COHERENT tt^ PRODUCTION 

In neutrino-nucleus interactions, there are two main mechanisms for production. The 
can result from the decay of a resonance, such as a A+ or A*^, that was produced in the 
primary interaction, or it can be produced coherently off of the entire nucleus. Coherent 
and resonant production have very different distributions for the pion angle with respect 
to the beam direction - the extra mass of the resonance tends to broaden the angular 
distribution, while the coherent pions pile-up in the forward direction. This fact can be 
used to fit the relative contribution of the two production mechanisms. 

The candidate events in the momentum reweighted MC are used to form three 
templates: one for resonant events, one for coherent events, and one for background. 
While the angular distribution for the coherent and resonant events are quite different, 
the coherent and background events are somewhat similar. So the templates are made 
in two dimensions: the first dimension is a function of angle (£';r(l — cos0)) and the 
second dimension is mass. The more complex angular function is used because it has 
a consistent shape for coherent events across all momenta at MiniBooNE energies. 
Variable binning is used such that the total number of MC events in each bin is approxi- 
mately equal. The number of bins in each projection is varied, independently, from 15 to 
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FIGURE 3. A comparison of raw and corrected Monte Carlo to data for various n kinematic distribu- 
tions: a) the opening angel between the two gammas, b) energy of the most energetic gamma, c) energy 
of the least energetic gamma, d) 77 mass. In all plots the corrected Monte Carlo is in as good or better 
agreement with the data than is the raw Monte Carlo. 



25, for a total of 121 fits, and the average fit parameters are used. Figure |4] shows the fit 
result plotted against data in the two projections. For the MiniBooNE flux and detector, 
with the NUANCE generator [4] providing the secondary interaction model, the fit finds 
that (19.5±1.1)% of all exclusive neutral current production is coherent. This should 
be compared to the raw MC which predicts 30% coherent for the MiniBooNE flux and 
detector. 

To determine the best overall production parameterization, the binned momentum and 
coherent fits are iterated. The iteration converges after only one round. In the oscillation 
analysis, the MC is corrected in both momentum and coherent fraction. 

In the coherent analysis, a number of possible sources of systematic error were inves- 
tigated including: choice of binning, background composition, momentum reweighting, 
neutrino flux, choice of analysis cuts, and detector model. Table [U lists the error esti- 
mation from each of these sources. By far, the largest source of systematic error is the 




detector model, which comes primarily from the uncertainty in the reconstructed energy 
scale. 

As a test of the model dependence of the coherent fraction, the data were refit after 
significant modifications of the MC model parameters. Five variations on the NUANCE 
model were explored: 

• The coherent axial mass assumed in the Rein-Sehgal coherent model ^ was 
decreased by a factor of three. This particular excursion was chosen because it 
is what would be needed to bring the predicted coherent cross section into 
agreement with the measured MiniBooNE coherent rate. While the normalization is 
given by the fit, this change does alter the distribution of kinematics for the coherent 
events, and therefore can change the fit result. 

• Diffractive events, which arise from coherent scattering off hydrogen targets, ac- 
count for 16% of all coherent events in MiniBooNE. In this variation diffractive 
events were removed from the coherent fit template. The diffractive contribution 



TABLE 1. Contributions to the systematic errors in the 
coherent fraction. Errors are given in percent coherent. 



Source Error (%) 



Binning 


0.21 


Background Model 


0.64 


Reweighting 


0.51 


Flux 


0.06 


Analysis Cuts 


0.51 


Detector Model 


2.34 


All Systematics 


2.54 



tends to be slightly less forward peaked than the coherent scatters off carbon, and 
hence their removal impacts the shape of the coherent fit template. 

• The axial mass, M^^, assumed in the Rein-Sehgal resonant model [6] was varied 
by ±25% of its default value. Altering the axial form factor parameter affects the 
resonant contribution at low Q^. 

• The nuclear model affecting resonant events was modified: 

- The binding energy in the Fermi Gas model was varied by 177% from 34 MeV 
to 60 MeV. 

- The Fermi momentum in the Fermi Gas model was varied by 72% from 
246 MeV/c to 423 MeV/c. 

The variations are quite large and specifically impact the resonant predictions at 
low Q^. By reducing the nuclear effects, one can test whether large changes to the 
Fermi Gas model, not from coherent scattering, can improve agreement to data in 
the most forward (or low Q^) region. 

• The A decay angular distribution was varied from the default model, which is from 
Rein and Sehgal [5]. The variations included isotropic decay in the center of mass 
frame, a pure spin 1/2 decay, and a pure spin 3/2 decay. The default model in 
NUANCE is isotropic decay and, in making the switch to the Rein and Sehgal 
decay distribution, a significant shift in the coherent fraction was observed. 

The fitted coherent fraction and fit confidence levels for each of these studies are given 
in Table |2l In addition, the fit confidence level for fits with the coherent fraction fixed to 
zero are also given. The zero coherent fits show clearly show that none of these variations 
prefer a coherent free production model as has been suggested by the K2K charged pion 
study [0]. 

TABLE 2. Average fitted coherent fractions and confidence levels for fits with several varia- 
tions of the cross section model, including the study of the angular distribution of the A decay. 
The default model's angular distribution is given by the model of Rein and Sehgal [5]. The 
confidence level is also given for fits where the coherent fraction is fixed to zero. 



Variation 


Coherent 




Avg. C.L. (%) 




Fraction (%) 


Coh. 


No Coh. 


Default Model 


19.5±L1 


5.97 


L8xlO-^*' 


Coherent 


19.1±1.1 


5.73 


3.8x10-'^ 


No Diffractive 


17.9±1.0 


11.63 


5.7x10"'^ 


M^" Resonant Hi 


17.9±L1 


3.27 


1.5x10-^"* 


M^" Resonant Lo 


2L1±L1 


5.00 


7.2x10-22 


Binding Energy Hi 


19.4±1.1 


5.85 




Binding Energy Lo 


19.6±1.1 


6.29 


2.2x10-^^ 


Fermi Momentum Hi 


18.2±1.1 


3.29 


1.3x10-15 


Fermi Momentum Lo 


2L0±L1 


4.24 


4.6x10-23 


Isotropic A Decay 


18.1±L1 


L88 


1.8x10-1*' 


Pure Spin 3/2 


20.8±1.0 


3.49 


1.7x10-1^ 


Pure Spin 1/2 


16.9±1.2 


0.01 


1.3x10-1'* 



CONCLUSIONS 



Neutral current Tt production is both a major potential background to the MiniBooNE 
oscillation analysis and an opportunity to make high impact measurements in neutrino 
cross sections (with the world's largest data set of 0.5 to 2 GeV neutrino interactions). 
The analysis described here has resulted in a direct measurement of production in 
the MiniBooNE detector, which is critically important for estimating the Tt^ misidenti- 
fication background to appearance. In addition, the coherent production, relative 
to total exclusive production, was measured and found to be (19.5 ±1.1 (stat) ±2.5 
(sys))%. 
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